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■ Abstract. After reviewing the definition of the heavy quark-antiquark potential in pNRQCD, we 

' discuss recent advances in the calculation. 
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<N : DEFINITION 



X 



The potential between a heavy quark and antiquark has been one of the first quantities 

to be studied in QCD: it is a privileged object for exploring the interplay of perturbative 

and non-perturbative QCD and the set in of confinement, and it plays a central role in 

^ I quarkonium physics [1]. Nowadays, the progress of perturbative and lattice calculations 

requires an accurate and rigorous definition of the potential in QCD, phenomenological 

and intuitive characterizations being no longer adequate. 

So, what is the QCD potential between a quark and antiquark with a large mass m? 

^ I One may first answer that the potential is the function V into the Schrodinger equation 

■ describing the quark-antiquark bound state 0: 

rn , 

O: E^=(^ + V](^, (1) 

O 

p being the momentum of the quark-antiquark pair in the centre-of-mass system and E 
its binding energy. Clearly, if Eq. ([T]) comes from a systematic expansion of QCD, it 
arises from at least a double expansion in p/mox rm (r being the inter-quark distance) 
I and in E r. Hence, rather than Eq. ([T]), we may expect that the QCD expansion would 

lead to 

E^=(pl + viyr) + ^-^ + ..\^, (2) 
\m m j 

where the . . . stand both for terms suppressed in the non-relativistic expansion m p/m 
or rm and for terms suppressed in Er, sometimes referred to as retardation effects (an 
example is the Lamb-shift). The above double expansion becomes an expansion in the 
heavy-quark velocity v once we note that in a non-relativistic system 1 /r ~ p ~ mv, and 
E ~ mv^, with v ^ 1 . 

How do the scales mv and mv^ originate in QCD? Let's consider the case of weakly- 
coupled bound states, i.e. states such that Aqcd is smaller than any of the scales m, mv 
or mv^. For these states we may use perturbation theory. Near threshold, the momenta 



as(l + as/v + ...) 

FIGURE 1. Resummed propagator near threshold. 



of the quarks are small compared to their masses, so that p/m ~ v <^ 1. Moreover, for 
certain sets of graphs, like those in Fig.[Tl the perturbative expansion breaks down when 
OCs ~ V. The summation of all (x^/v contributions leads to the appearance of a bound-state 
pole of order mv^ ~ mo^ in the resummed propagator. 

These scales get entangled in a typical amplitude. An example is provided by the 
annihilation diagram of Fig. [2l Assuming that the incoming quarks are near threshold, 
the different gluons entering the diagram are characterized by different scales. The 
annihilation gluons have a typical energy of order m, sometimes also called "hard scale"; 
binding gluons, also called "soft", have the momentum of the incoming quarks, which 
is of order mv, and "ultrasoft" gluons, sensitive to the intermediate bound state, have 
energies of the order of the binding energy, i.e. mv^. 
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FIGURE 2. Annihilation diagram contributing to the quarkonium decay width. 

In order to disentangle the different scales, it is convenient to enforce an expansion 
in the ratios of low-energy scales over large-energy scales at the Lagrangian level; this 
corresponds to substituting QCD with low-energy Effective Field Theories (EFTs) 10]. 
The ultimate EFT that follows from QCD by integrating out all energy scales but mv^ is 
potential NRQCD (pNRQCD). The general form of the Lagrangian density of pNRQCD 
is 

if = [ d^rS^ ( ido- — - (r, {x) - ^^^'^^ + . . . V + ultrasoft contributions, 
J \ m ml 

(3) 

where S stands for a color-singlet quarkonium field, ji is the cut-off of the EFT, and 
"ultrasoft contributions" include all degrees of freedom which are ultrasoft (they may be 



gluons, or light quarks or other degrees of freedom). The ukrasoft contributions are 
typically suppressed with respect to the part of the pNRQCD Lagrangian displayed 
in Eq. (|3]). Hence, the equation of motion of the color-singlet quarkonium field is 
exactly Eq. Q and we may identify V^^\r,ii)+V^^\r,ii)/m + . . . with the heavy-quark 
potential. 

In summary, EFTs provide the following definition of the potential: the potential is 
a Wilson coefficient of the EFT obtained by integrating out all degrees of freedom 
but the ultrasoft ones, it undergoes renormalization, develops a scale dependence and 
satisfies renormalization group equations, which eventually allow to resum potentially 
large logarithms. 



THE PERTURBATIVE POTENTIAL 



We consider the static potential Vs^\ This is obtained by integrating out soft gluons 
from static QCD. Soft gluons are those associated with the scale 1/r. At short distances, 
1/r 3> Aqcd, soft gluons may be calculated in perturbation theory. If also the ultra- 
soft scale, i.e. the potential itself, is larger than Aqcd? then, besides the color-singlet 
quarkonium field, ultrasoft degrees of freedom include ultrasoft gluons and the color- 
octet quarkonium field. The matching leading QCD to pNRQCD may be done in per- 
turbation theory, see Fig. [31 Sometimes it may be useful to choose the QCD Green's 
function in a gauge invariant fashion. A popular choice is the static Wilson loop. 
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FIGURE 3. Matching condition for pNRQCD. On the left-hand side the four-fermion Green's function 
in QCD, on the right-hand side the singlet propagator and the first ultrasoft correction in pNRQCD. The 
single continuous line stands for a singlet propagator, the double line for an octet propagator and the curly 
line for a chromoelectric correlator The coupling of the gauge fields with the quarkonium (circle with a 
cross) is a chromoelectric dipole vertex. 

The matching fixes the potential and the other Wilson coefficients of the EFT. The 
matching condition for the singlet static potential reads 



lim — ln< 



Nr 



(Tr(E(0 ■E(0)))(Ai) + . . ., (4) 



where the box stands for the static Wilson loop of dimension rxT, vj^^ for the static 
octet potential, Va for the electric-dipole matching coefficient, E for the chromoelectric 



field and Nc — 3 for the number of colors. The left-hand side of Eq. (HI) is known at two 
loops [11,3, [!,[§]. At three loops the static Wilson loop contains a term proportional to 
cc^/r X In tts, which has been calculated in jvilil]. 

In order to determine the matching coefficients vj^^ and Va that enter in Eq. dU besides 

, we need two further matching conditions. The static octet potential Vo^^ has been 
calculated up to two loops by matching it to a static Wilson loop with color matrices in 
the initial and final states [|9i]. This gives rise to a matching condition similar to Eq. dH) 
i]: 

limi;ln<Q> = vi°)(r,Ai) + .... (5) 
The matching for Va is described in Fig. H it gives 

VA{r,^) = l + 0{a^). (6) 



+ 




FIGURE 4. Matching condition for Va- 

The last ingredient needed in order to calculate Eq. dH) is the chromoelectric cor- 
relator (Tr(E(?) -£(0))), where Wilson lines connecting the chromoelectric fields are 
understood. This has been calculated at order in jioll . 

Since the static Wilson loop is fully known at two loops, the matching condition 
dH) provides the static singlet potential at two loops. Moreover, since the static Wilson 
loop is independent of /i, the right hand-side of Eq. @) should also be /^-independent. 
Therefore, the logarithmic dependence of the static potential may be extracted by not- 
ing that the Inr/i, In^rji, ... terms in Vs^^ have to cancel against the ln(yj'^'' — vj-^'')/ ji, 

ln2(vi°) - vJ^Vm, - Inr^i, In^ r^i, ... terms in Tj^e-''^^''"'-^^'"') (Tr(E(0 -£(0))). This 

Jo 

leads to a great simplification in the calculation of the logarithmic dependence of the 
static potential: the logarithmic contribution at N-^LO and the single logarithmic contri- 
bution at N^LO may be extracted respectively from a one-loop and two-loop calculation 
in the EFT. Finally, we note that the solutions of the renormalization group equations 
allow the calculation and resummation of all logarithmic contributions of a given type 
(e.g. leading logarithms of the type al x (a^ In /ir)'\ next-to-leading logarithms of the 
type x ( In /ir)" and so on). 



The presently most accurately known fixed-order expression of the static singlet 
potential is 
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.2^3, 



(7) 



where Cp = Tf{N^ - l)/Nc, Ca = Nc, Tp = 1/2, /3o = 1 1Ca/3 -4r^R//3, /3i = 34Cl/3 
—2QCATpnf/?i — ACpTFnf, nf is the number of (massless) flavors, Je is the Euler 
constant and Ok is the strong coupling constant in the MS scheme. The coefficients a\, 
a2, and a\ stand for 
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(9) 
(10) 



«i + 2r£/3o + n/ f - — + - In 2 j 
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(11) 



The coefficient a\ was calculated in [11], the coefficient a? in [jsj 3, HI |6l], the term 
proportional to a^/r x In r/i in [SB], the coefficient in [l2,[l3l] and the coefficient 
«4 in [I13I1 . The coefficients 03 and 04 are only partially known (see [|13n for discussion 
and references). The leading logarithmic contributions have been resummed to all orders 
infl?]. 

Expression ([7]) shows explicitly the Wilson coefficient nature of the static potential. 
It shows a scale dependence, which comes from the renormalization, and it satisfies 
renormalization group equations, which allow to resum potentially large In r/i terms. 
Also large contributions of the renormalon type may be analyzed in the EFT framework. 



By summing Eq. (Q to the ultrasoft contributions we get back the static Wilson loop, 
i.e. the energy between two static sources in QCD. This reads 
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«2+(y + 4T^)/3o' + 2r£(2ai/3o + /3i; 



16;r2 Q^(l/r) ■ 
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af In^ ^ii^^ +4ln^^i^iilM +a- 



(12) 



This quantity may be compared with the short-distance behaviour of the static Wilson 
loop provided by lattice calculations, see for instance [il4i [l5ll. 



THE NON-PERTURBATIVE POTENTIAL 



At large distances, 1/r ~ Aqcd, due to confinement, ultrasoft effective degrees of 
freedom may only be colorless objects. If Goldstone bosons are neglected, the color- 
singlet quarkonium field S turns out to be the only dynamical degree of freedom at 
scales lower than Aqcd fs*]. The static singlet potential is then simply given by 



(0) 



lim — ln< 



(13) 



A recent lattice determination is shown in Fig. [51 




r [fm] 

FIGURE 5. Lattice determination of the right-hand side of Eq. ( fT3T l. from llf 

Recently, and for the first time, the leading relativistic correction to the static potential 
has been calculated on the lattice. The existence of a possibly large non-perturbative \/m 

potential, Vs'^\ was first pointed out in [ITtIi . v/^-* may be written as a static Wilson loop 
with two chromoelectric field insertions on the same quark line: 



(1) 



m 
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= I dtti\ 

2m Jo I 1 



>■ 



(14) 



The corresponding lattice determination is shown in Fig. [61 
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FIGURE 6. Lattice determination of the right-hand side of Eq. (O, from IH] . 

Note that, in accordance to power counting arguments, in the long-range, the 1/m 
potential may be as large as the static potential and contribute with it to the leading- 
order potential [vf]. 

By the same collaboration, spin-dependent 1 / potentials have been calculated on 
the lattice with unprecedented precision. Expressions for the spin-dependent potentials 
in terms of static Wilson loops and field- strength insertions have been derived in [1^ 20ll . 
These have been used to obtain the lattice results shown in Fig.|71 
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FIGURE 7. Lattice determination of the spin-dependent 1 /rrt^ potentials, from 1 16]. The potentials Vj' 
and V2 are spin-orbit potentials, the potential V3 is the tensor potential and the potential V4 the spin-spin 
potential. 

In the long range, the spin-orbit potentials show, for the first time, deviations from 
the flux-tube picture of chromoelectric confinement. Since a fully consistent renormal- 



ization of the EFT operators is still missing in the lattice analysis, it may be premature 
to draw any definitive conclusion. However, progress has been made recently in this di- 
rection. In [21], the non-perturbative renormalization of the chromomagnetic operator 
in the Heavy Quark Effective Theory, which crucially enters in all spin-dependent po- 
tentials, has been performed for the first time. A proper operator renormalization is also 
crucial in order to verifvan exact relation among the spin-dependent potentials required 
by Lorentz invariance [|22l.l23ll. which was checked in lfl6ll at the few percent level. 



CONCLUSIONS 

Non-relativistic EFTs provide a rigorous definition of the potential between a heavy 
quark and antiquark (see [24] for systems made by two or three heavy quarks). In the 
perturbative regime, the potential is a key ingredient for precision calculations of several 
threshold observables. In the non-perturbative regime, it can be calculated on the lattice; 
the corresponding EFT, pNRQCD, may provide lattice studies with an alternative to 
more traditional EFTs with heavy quarks, like NRQCD. 
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